Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/M  R/6790-  -07-9043 


Design  of  a  Compact,  Optically 
Guided,  Pinched,  Megawatt  Class 
Free-Electron  Laser 


Phillip  Sprangle 
Joseph  PeSano 

Beam  Physics  Branch 
Plasma  Physics  Division 


Bahman  Hafizi 

Icarus  Research,  Inc. 
Bethesda,  Maryland 


June  8,  2007 

2007091 2340 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


PuWie  reporting  human  lor  this  collection  of  infomaiion  is  estimated  to  average  1  hour  per  response,  including  me  lime  (of  reviewing  instructions,  searching  existing  data  sources,  gathonng  ana 

maintaining  the  data  needed,  and  compieiing  and  reviewing  ttiia  collection  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  coHection  of  inlonmabon.  Including 
suggestions  lor  reducing  this  burden  to  Department  oJ  Delense,  Washington  Headquarters  Seivices.  Directorate  lor  information  Operations  and  Reports  (0704-0168),  1215  Jefferson  Davis  Highway, 
Suite  1204.  Arlington,  VA  22202*4302.  Respondents  should  he  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  fading  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number  PLEASE  DO  MOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS, 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

08-06-2007  Final 

3.  DATES  COVERED  (From  -  To) 

May  2004  -  May  2006 

4,  TITLE  AND  SUBTITLE 

Design  of  a  Compact,  Optically  Guided,  Pinched, 

Megawatt  Class  Free- Elec  iron  Laser 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Phillip  Sprangle,  Joseph  Pefiano.  and  Bahman  Hafizi* 

5d.  PROJECT  NUMBER 

67-8915-07 

5e,  TASK  NUMBER 

5f,  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

Naval  Research  Laboratory,  Code  6790 

4555  Overlook  Avenue.  SW  NRL/MR/6790-07-9043 

Washington.  DC  20375-5320 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSOR  /  MONITOR'S  ACRONYM(S) 

Office  of  Naval  Research 
One  Liberty  Center 
875  North  Randolph  St. 

Arlington,  VA  22203-1995 

12.  DISTRIBUTION  f  AVAILABILITY  STATEMENT 


QNR 


11  *  SPONSOR  /  MONITOR'S  REPORT 
NUMBER(S) 


Approved  for  public  release:  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 

♦Icarus  Research,  Inc.,  P.O,  Box  30780,  Bcthesda,  MD  20874-0780 


14,  ABSTRACT 

A  conceptual  design  for  o  compact  megawatt  class  FEL  operating  within  atmospheric  transmission  windows  is  presented.  The  proposed  FEL 
consists  of  an  optically  guided,  pinched  amplifier  configuration  driven  by  an  RF  linac.  The  gain  length,  efficiency,  electron  pulse  slippage  and  the 
distance  between  the  wiggler  and  first  relay  mirror  are  determined  for  a  megawatt  class  design.  Of  particular  concern  in  the  design  is  the  overall 
length  of  the  optical  system,  Le.t  wiggler  length  and  distance  to  the  first  relay  mirror  En  the  present  design  the  wiggler  length  is  -2  meters  and 
the  distance  between  the  first  relay  mirror  and  the  wiggler  is  determined  by  the  average  intensity  damage  threshold  on  the  mirror.  By  focusing 
the  electron  beam,  the  optical  beam  can  be  pinched  upon  exiting  the  wiggler.  The  pinched  optical  beam  has  a  reduced  Rayleigh  length  which 
permits  the  first  relay  mirror  to  be  relatively  close  to  the  wiggler.  By  pinching  the  optical  beam  and  employing  grazing  incidence  the  first  relay 
mirror  can  be  located  with  ~3  meters  of  the  wiggler.  It  is  shown  that  frequency  detuning  can  more  than  double  the  FEL  efficiency.  In  addition, 
electron  pulse  slippage  is  shown  to  be  substantially  reduced  in  a  high-gain  amplifier. 


15.  SUBJECT  TERMS 


Megawatt  FEL 
Optically  guided 

High-gain 

Compact 

Pinched 

Amplifier 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Phillip  Sprangle 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

UL 

30 

19b.  TELEPHONE  NUMBER  (tnefode  area 
code) 

(202)  767-3493 

l 


Standard  Form  298  (Rev.  8-98) 

Pr*ftcHb*d  by  ANSI  Sid,  Z39.1S 


Final  Technical  Report 


Design  of  a  Compact,  Optically-Guided,  Pinched, 
Megawatt-Class  Free-Electron  Laser 


Prepared  for 

HIGH-ENERGY  LASER  JOINT  TECHNOLOGY  OFFICE 


For  the  period 
May  2004  to  May  2006 


Prepared  by 

Phillip  Sprangle  (Principal  Investigator) 

Joseph  Pehano 
Bahman  Hafizi 

Icarus  Research,  Inc.,  P.O.  Box  30780.  Bethesda.  MD  20824-0780 


Manuscript  approved  May  1 4.  2007 


CONTENTS 


SUMMARY . . iv 

I.  Introduction . I 

II.  Laser  Propagation  in  a  Maritime  Environment . .  2 


i)  Maritime  Atmosphere  Model 

ii)  HEL  Propagation  Simulations 

III.  Optically  Guided  FEL  Amplifier . . 5 

i)  Operating  Wavelength 

ii)  Gain  Length 

iii)  FEL  Efficiency 

iv)  Optical  Guiding 

v)  Saturated  Peak  and  Average  Optical  Power 

vi)  Electron  Pulse  Slippage  (Lethargy) 

vii)  Pinched  Optical  Beam 

viii)  Optical  Beam  Quality  and  Beam  Spreading 

ix)  Distance  Between  First  Relay  Mirror  and  W'iggler 

x)  Frequency  Detuning  and  Efficiency  Enhancement 


IV.  Discussion  and  Conclusions . 10 

Acknowledgments . 12 

Appendix . 13 

REFERENCES . 15 

FIGURES . 18 


iii 


TABLES 


24 


SUMMARY 

A  conceptual  design  for  a  compact  megawatt  class  FEL  operating  within  atmospheric 
transmission  windows  is  presented.  The  proposed  FEL  consists  of  an  optically  guided, 
pinched  amplifier  configuration  driven  by  an  RF  linac.  The  gain  length,  efficiency, 
electron  pulse  slippage  and  the  distance  between  the  wiggler  and  first  relay  mirror  are 
determined  for  a  megawatt  class  design.  Of  particular  concern  in  the  design  is  the  overall 
length  of  the  optical  system,  i.e.,  wiggler  length  and  distance  to  the  first  relay  mirror.  In 
the  present  design  the  wiggler  length  is  -  2  meters  and  the  distance  between  the  first  relay 
mirror  and  the  wiggler  is  determined  by  the  average  intensity  damage  threshold  on  the 
mirror.  By  focusing  the  electron  beam,  the  optical  beam  can  be  pinched  upon  exiting  the 
wiggler.  The  pinched  optical  beam  has  a  reduced  Rayleigh  length  which  permits  the  first 
relay  mirror  to  be  relatively  close  to  the  wiggler.  By  pinching  the  optical  beam  and 
employing  grazing  incidence  the  first  relay  mirror  can  be  located  within  -  3  meters  of  the 
wiggler.  It  is  shown  that  frequency  detuning  can  more  than  double  the  FEL  efficiency. 
In  addition,  electron  pulse  slippage  is  shown  to  be  substantially  reduced  in  a  high  gain 
amplifier. 
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1.  Introduction 

The  free  electron  laser  (FEL)  is  capable  of  producing  high  average  power  at  high 
efficiency  without  the  conventional  thermal  management  and  waste  issues  associated 
with  other  laser  systems  [1  ].  In  addition,  the  operating  wavelength  can  be  chosen  for 
optimized  propagation  in  a  maritime  environment  [2],  These  unique  features  make  the 
FEL  a  leading  candidate  for  naval  directed  energy  applications  [3-6], 

Propagation  of  high  power  laser  beams  in  the  atmosphere  involves  a  complex, 
linear  and  nonlinear  interaction  between  various  physical  processes.  The 
theoretical/numerical  model  used  in  this  study  includes  the  effects  of  aerosol  and 
molecular  scattering,  aerosol  heating  and  vaporization,  thermal  blooming  due  to  aerosol 
and  molecular  absorption,  atmospheric  turbulence,  and  beam  quality.  These  processes  are 
modeled  in  a  fully  three-dimensional  and  time-dependent  manner.  It  is  found  that 
aerosols,  which  consist  of  water,  sea  salt,  organic  matter,  dust,  soot,  biomass  smoke, 
urban  pollutants,  etc.,  are  particularly  important  because  they  result  in  laser  scattering, 
absorption  and  enhanced  thermal  blooming. 

In  this  report  we  present  parameters  for  a  conceptual  design  of  a  compact, 
megawatt-class,  optically  guided,  free-electron  laser  amplifier.  There  are  a  number 
important  characteristics  and  advantages  of  the  proposed  configuration.  In  an  optically- 
guided  FEL  amplifier  the  filling  factor,  i.e,.  the  ratio  of  the  electron  beam  to  the  optical 
beam  cross-sectional  area,  remains  relatively  constant  throughout  the  wiggler  region,  and 
the  radiation  amplitude  increases  exponentially  until  saturation.  The  optical  beam  upon 
exiting  the  wiggler  can  be  pinched  by  focusing  the  electron  beam  with  external  focusing 
fields.  The  Rayleigh  range  of  the  pinched  optical  beam  can  be  short.  As  a  result  the  first 
relay  mirror  can  be  placed  close  to  the  wiggler  without  exceeding  the  intensity  damage 
threshold.  The  first  relay  mirror  is  a  grazing  incidence  reflector.  The  FEL  efficiency  can 
be  increased  by  tapering  the  wiggler.  Alternatively,  in  a  uniform  wiggler,  enhanced 
efficiency  can  be  obtained  by  frequency  detuning  the  input  signal.  For  the  power  levels 
of  interest  we  consider  FEL  operation  at  wavelengths  1  pm  and  2.141  pm,  i.e.,  two 
transmission  windows  associated  with  atmospheric  water  vapor.  Because  of  propagation 
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and  generation  advantages  as  well  as  for  eye  safety  reasons,  operation  at  2.141  pm  may 
be  more  advantageous  [4]. 

Recently  a  high-gain  FEL  amplifier  operating  in  the  optically  guided  regime  has 
been  experimentally  demonstrated  at  the  Brookhaven  National  Laboratory  [8].  Our 
analysis  of  this  experiment  is  in  close  agreement  with  the  experimental  results  and  will  be 
discussed  in  more  detail  in  Sec.  IV. 

II.  Laser  Propagation  in  a  Maritime  Environment 

The  optimum  wavelength  for  efficient  HEL  propagation  depends  on  the 
atmospheric  conditions  and  a  number  of  inter-related  physical  processes  which  include; 
thermal  blooming  due  to  aerosol  and  molecular  absorption  [9],  turbulence  [10],  aerosol 
and  molecular  scattering  [11],  thermal  scattering  due  to  heated  aerosols,  and  aerosol 
heating  and  vaporization  [12-15].  The  relat  ive  importance  of  these  processes  depends  on 
the  parameters  of  the  atmospheric  environment  which  can  vary  significantly  depending 
on  location  and  time. 

Atmospheric  environments  contain  various  types  and  concentrations  of  aerosols. 
Aerosols  can  absorb  laser  energy  and  conductive ly  heat  the  surrounding  air,  resulting  in 
an  increase  in  thermal  blooming  of  the  HEL  beam  [2,16].  In  general,  aerosols  consist  of 
hygroscopic  and  non-hygroscopic  particles  of  various  sizes  and  chemical  compositions. 
Oceanic  aerosols  consist  of  sea  salt,  water,  and  organic  material.  Non-hygroscopic 
aerosols  are  composed  of  dust,  soot,  biomass  smoke,  and  other  carbon-based  compounds. 
These  aerosols  typically  have  much  larger  absorption  coefficients  than  water-based 
aerosols.  While  they  are  normally  present  in  continental,  rural  and  urban  environments, 
dust  aerosols  can  also  be  present  in  maritime  environments  hundreds  of  miles  from  shore 

[17]. 

i)  Maritime  Atmosphere  Model 

The  Advanced  Navy  Aerosol  Model  (ANAM)  is  used  to  model  the  near-surface 
maritime  environment  [18],  The  ANAM  aerosol  distribution  is  comprised  of  various 
modes  which  represent  aerosols  of  different  compositions  and  sizes.  These  aerosol 
modes  will  absorb  laser  energy  and  vaporize  at  different  rates.  The  aerosol  size 
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distribution  function,  F(R.t)  =  ^  F,(Rj ) .  where  R  is  the  aerosol  radius,  is  represented 

j  =  o 

as  a  superposition  of  five  “modes''  with  each  mode  representing  aerosols  with  a  particular 
physical  composition  and  origin.  The  aerosol  distribution  can  evolve  with  time,  due  to 
vaporization  for  example,  the  aerosol  absorption  and  scattering  coefficients  are  also  time- 
dependent.  Mode  0  represents  dust  particles  of  continental  origin,  mode  1  represents 
water-soluble  aerosols,  and  modes  2-4  represent  marine  aerosols  (sea  salt  and  water)  that 
result  from  different  processes. 

We  use  MODTRAN  to  calculate  the  molecular  absorption  and  scattering 
coefficients  [19].  Figure  1  plots  the  gross  scattering  and  absorption  coefficients  due  to 
aerosol  and  molecular  contributions  in  the  water  vapor  transmission  windows, 

X  =  1 .045, 1 .625  and  2.141  pm .  The  gross  scattering  and  absorption  coefficients  that  we 
obtain  are  comparable  with  in  situ  measurements  [20. 21  j. 

ii)  HEL  Propagation  Simulations 

To  simulate  HEL  propagation  we  use  the  High  Energy  Laser  Code  for 
Atmospheric  Propagation.  HELCAP  [2,22].  developed  at  the  Naval  Research  Laboratory. 
HELCAP  includes  the  effects  of  i)  aerosol  and  molecular  scattering,  ii)  aerosol  heating 
and  vaporization.  Hi)  thermal  blooming  due  to  both  aerosol  and  molecular  absorption,  iv) 
atmospheric  turbulence,  and  v)  laser  beam  quality,  i.e..  all  of  the  effects  discussed  in  the 
preceeding  section.  It  is  the  first  HEL  propagation  model  which  integrates  all  these 
physical  processes  in  a  fully  three-dimensional,  time-dependent  manner.  In  modeling  the 
aerosol  effects,  w:e  account  for  the  aerosol  distribution  and  the  various  aerosol  modes 
(water-based,  dust.  soot.  etc.). 

In  HELCAP  the  laser  electric  field  is  represented  as 
E  =  ^(x,>’,z,f)exp[  i(coz/c  -  (ot)\  eT  /  2  +  c.c. ,  where  a)  =  2nd  X  is  the  laser 

frequency.  eT  is  a  unit  polarization  vector  in  the  *  direction.  A(x,y.z,t )  is  the  complex 

laser  amplitude  and  the  laser  intensity  is  !  =  cAA‘ /8 n .  HELCAP  solves  a  nonlinear 
Schrddinger-like  equation  which  has  the  form 
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d  A 
dz 


=  ~W\A 
2m  x 


+  dn  JH ) 


I** + » 


A , 


where  a  =  am  +  aA  is  the  total  absorption  coefficient,  /?  =  /?„,+  is  the  total 
scattering  coefficient,  and  Snr  and  dnTB  denote  the  refractive  index  variation  due  to 
atmospheric  turbulence  and  thermal  blooming  respectively.  am  ( a A  )  is  the  molecular 
(aerosol)  absorption  coefficient.  ( fiA  )  is  the  molecular  (aerosol)  scattering 
coefficient.  The  quantities dw, ,  dnTH ,  a ,  and  f3  are  space  and  time-dependent  and 

determined  self-consistently  in  the  presence  of  aerosol  heating  and  vaporization. 

The  transmitted  power  at  the  source  is  denoted  by  P  and  the  laser  is  focused  onto 
a  remote  target  at  a  range  of  5  km  (see  Fig.  2).  The  target  is  taken  to  be  circular  with  an 
area  of  100  cm2.  The  propagation  direction  is  along  the  z-axis  and  a  uniform  transverse 
wind,  with  velocity  Vw  =  5  m/sec  is  directed  along  the  y-axis.  Atmospheric  turbulence  is 

modeled  by  a  Kolmogorov  spectrum  with  structure  constant  C;  =  I0"15  m  2 ' .  The 

pointing  jitter  associated  with  the  laser  beam  is  taken  to  have  an  angular  spread  of  2  prad 
and  a  white  noise  temporal  spectrum. 

Figure  3  plots  the  average  power  on  target  versus  the  transmitted  power,  P ,  for 
the  three  wavelengths  of  interest.  Our  results  show  that  for  a  maritime  environment,  the 
optimum  wavelength  depends  on  the  transmitted  power.  For  P  <  1.5  MW,  propagation 

is  mostly  affected  by  aerosol  scattering  and  the  average  power  on  target  increases  with 
P  .  In  this  regime,  the  1 .625  p  m,  and  2.141  pm  wavelengths  provide  slightly  greater 
power  on  target  than  1 .045  p  m.  This  is  due  to  the  lower  aerosol  scattering  coefficient 
associated  with  the  longer  wavelengths.  For  P  <  1.5  MW,  the  propagation  efficiency  is 
roughly  50%  for  the  three  wavelengths  considered.  For  example,  P  =  1  MW  results  in 

(Png#)  -  0.55  MW  for  X  =  1.625  pm  and  X  =  2.141pm.  However,  for  P  >1.5  MW, 

thermal  blooming  becomes  important.  In  this  high  pow:er  regime  the  optimum 
wavelength  is  1 .045  pm  due  to  the  lower  molecular  absorption  coefficient  in  that  water 

vapor  window.  For  P  =  3  MW.  (fjarget)  ~  1  MW  for  1.045  pm  while  -  0.8  MW 

for  1.625  pm  and  2.141  pm. 
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For  the  power  levels  of  interest.  P  <  1 .5  MW  average  power,  it  may  be 
advantageous  to  operate  the  FEL  at  2.141  pm  .  For  maritime  propagation  there  is  no 
distinct  advantage  in  operating  the  FEL  at  either  1 .045  pm  or  1 .625  pm.  However,  for 
P  >1.5  MW ,  operating  at  1  pm  leads  to  greater  propagation  efficiency. 

III.  Optically  Guided  FEL  Amplifier 

Using  the  Source  Dependent  Expansion  (SDE)  formulation  [23],  expressions  for 
the  optical  growth  rate,  wavenumber  shift,  spot  size  and  optical  curvature  of  the  optical 
beam  in  an  FEL  amplifier  are  obtained  in  the  Appendix.  These  results  are  used  to  design 
a  compact  MW  FEL  in  the  optically  guided,  high  gain  regime,  i.e.,  having  a  matched 
optical  spot  size.  In  addition,  focusing  the  electron  beam  at  the  wiggler  exit  is  shown  to 
result  in  a  pinched  optical  beam  with  a  shortened  Rayleigh  range.  A  schematic  diagram 
of  the  FEL  amplifier  configuration  is  shown  in  Fig.  4.  Finally,  it  is  shown  that  by 
frequency  detuning  the  FEL  efficiency  can  be  more  than  doubled.  The  following 
discussion  motivates  the  design  parameters  listed  in  Tables  I  -  V.  For  illustrative 
purposes,  the  numbers  given  in  the  text  correspond  to  a  resonant  wavelength  of 
X  =  2.141  fjm  with  Am  =  0  unless  otherwise  stated. 

i)  Operating  Wavelength 

In  the  present  design  the  resonant  operating  wavelength  of  a  linearly  polarized 
optical  beam  is  X  =  Xu,  (1  +  K1  /2)/2y 2 .  where  y  =  fT^MeVJ/0.5 11  +  1  is  the  total 
relativistic  mass  factor.  Design  parameters  for  X  -  1  pm  and  X  =  2.141  /tm  are  listed  in 
Tables  I-V. 

ii)  Gain  Length 

The  power  gain  length  in  the  optically  guided  regime,  in  the  absence  of  frequency 
detuning  (Aw  =  0).  is 


V' 


=  14cm,  (1) 
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where  /  =  R';  / Rj  <1  is  the  filling  factor,  v  =  /  JkA]/17  is  Budker’s  parameter,  Ih  is 
the  electron  beam  current  in  kA,  K  =  q  BiV  12k  me'  =  1 .9  is  the  wiggler  strength 
parameter.  Bw  =  5kG  is  the  wiggler  field  and  =  4cm  is  the  wiggler  period.  The 
transverse  profiles  of  the  optical  and  electron  beam  are  Gaussian  with  radii  RL  and  Rh , 
respectively  [see  Eqs.  (A2)  and  (A4)]. 


iii)  FEL  Efficiency 

In  the  optically  guided,  high  gain  regime  the  efficiency  of  conversion  from 
electron  beam  power  to  optical  power  is  given  by  [3,24] 


n  = 


A.M  _  /,',(2  +  3/)',J 

V 

V2 

K 

2a  '  (1  +  2  f)2 

j  j 

+  K2/2 

=  0.8%, 


(2) 


Here,  A k  is  the  wavenumber  shift,  defined  in  Eq.  (A3).  Free-electron  laser  efficiency 
can  be  increased  by  tapering  the  wiggler  [5,  24,  25]  or  by  frequency  detuning  [3,24]. 


iv)  Optical  Guiding 

The  FEE  optical  beam  can  undergo  refractive  guiding  within  the  wiggler.  An 
envelope  equation  for  the  optical  beam  spot  size  can  be  obtained  by  using  the  results  in 
the  Appendix,  and  the  condition  for  a  matched  optical  beam  can  be  obtained.  For  optical 
guiding  the  power  gain  length  is  nearly  equal  to  the  free  space  Rayleigh  range,  that  is 

(2  +  3/)l/2  „ 


L,  = 


2fU2(l  +  2f) 


*  R  * 


(3) 


where  ZR  =  kR)J  k  =  22cm  is  the  Rayleigh  range.  The  optical  guiding  condition  in  Eq. 
(3)  can  be  rewritten  in  terms  of  the  beam  current. 


!h[kA\  =  1.1x10* 


2  f’(  1  +  2/)6 

4  (l  +  K2 12) 

<2  +  3  ff 

K  J 

y2  K2 

=  lkA . 


(4) 


In  addition,  the  centroid  of  the  FEL  optical  beam  can  be  steered  by  gradually  bending  the 
electron  beam. 
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v)  Saturated  Peak  and  Average  Optical  Power 

The  peak  optical  power  at  saturation  is  P  =  rj  I hVh  =1  GW  ,  where  Vh  =81  MV 
is  the  electron  beam  voltage.  The  average  power  at  saturation  is 

(P)  =  nDIh  Vb  =  1.5 MW .  (5) 

where  D  =  1 .5  x  1 0‘3  is  the  duty  factor  and  //  =  1 .2  %  ( A  co!  on  =  - 1 .5  %  ). 


vi)  Electron  Pulse  Slippage  (Lethargy) 

In  the  high-gain  regime  electron  slippage  (lethargy)  is  shown  to  be  significantly 
less  than  in  the  low  gain  regime  and  is  not  a  concern  for  this  design.  The  slippage  length, 
i.e.,  the  separation  between  the  optical  and  electron  pulse,  is  S  =  (vK  -  v_.)Z,„  /c,  where 

\  K  =  dcoldk  is  the  optical  pulse  group  velocity,  v.  =  c(l  -  \!2y:)  is  the  axial  electron 

beam  velocity,  and  y.  =  y/fi  +  12) .  Making  use  of  the  high  gain  FEL  dispersion 

relation  [3,24],  the  group  velocity  is  v  =  c(l  -  1/3 y)  -  {?J x Rt  )2  / 2),  and  the  slippage 
length  is 


S  =  AN.. 


[1 

2\ 

l1 ' 

[xRlJ 

(6) 


where  Nw  -  Lv  !  Av  is  the  number  of  wiggler  periods.  The  well  known  slippage  length 

in  the  low  gain  regime  is  A  Nw  .  Slippage  in  the  high  gain  regime  is  reduced  by  more 
than  a  factor  of  three. 


vii)  Pinched  Optical  Beam 

The  intensity  on  the  first  relay  mirror  must  be  below  a  damage  threshold  level. 

To  further  reduce  the  intensity  on  the  mirror  the  optical  beam  can  be  pinched  at  the 
wiggler  exit  in  order  to  shorten  the  Rayleigh  range  [3-5].  By  employing  external 
foeusing  fields  near  the  wiggler  exit  the  electron  beam  and  in  turn  the  optical  beam  can 
be  pinched.  Making  use  of  grazing  incidence  geometry  the  intensity  on  the  mirror  can  be 
further  reduced. 
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The  spot  size  of  the  optical  beam  R,(z)  [obtained  by  solving  Eqs.  (A5),  (A6)  and 
(A8)]  and  the  radius  of  the  electron  beam  Rb(z)  are  plotted  as  functions  of  the 

propagation  distance  z  in  Fig.  5.  The  entrance  to  the  wiggler  is  at  z  =  0  and  extends 
toz  =  Lw  =  1 .5  m  .  In  Fig.  5(a)  the  radius  of  the  electron  beam  is  constant  and  the  optical 

beam  is  guided  through  the  wiggler  with  a  constant  spot  size.  Upon  exiting  the  wiggler, 
the  optical  beam  diffracts.  In  the  example  of  Fig,  5(b)  the  electron  beam  is  focused  by 
external  fields  applied  near  the  wiggler  exit.  This  pinches  the  optical  beam  and  as  a 
result,  the  optical  beam  diffracts  more  rapidly  than  in  the  unpinched  example. 

viii)  Optical  Beam  Quality  and  Beam  Spreading 

It  is  common  practice  to  characterize  the  higher  order  modal  content  of  a  laser 
beam  by  a  beam  quality  parameter  denoted  by  M2 .  The  quantity  M2  >  1  is  a  “times 
diffraction-limited”  parameter  which,  for  a  fundamental  Gaussian  beam,  is  unity.  Upon 
exiting  the  FEL  interaction  region,  i.e.,  wiggler.  the  spatial  spreading  of  the  optical  beam 
is  given  by 

W\z)  =  ~\r2I{r,z)2nrdr  =  W2  +  M4  02D(z  -  LJ2 ,  (7) 

*  0 

where  fV  is  the  laser  beam  spot  size,  P  is  the  power  and  I  is  the  intensity.  In  Eq.(  7),  Wn 
is  the  minimum  spot  size  at  z  =  Lw ,  0D  =  X i{izW0)  is  the  diffraction  (spreading)  angle 
of  an  ideal  Gaussian  beam  and  M1  is  a  measure  of  the  beam  quality.  For  a  non-ideal 

optical  beam  the  diffraction  angle  is  M2  8n  at  the  exit  of  the  wiggler.  Using  the 
paraxial  wave  equation  in  Sec.  II,  it  can  be  shown  that  the  general  expression  for  the 
evolution  of  M 2  is 

M2{z)  = 

(8) 


4  x2 


\  '35 

jlr'dr  || 


a 


0  V 


dr 


ydrj  j 


■dr  -  |  \l^r2dr 


i  ■  l 
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where  the  complex  laser  amplitude  has  been  expressed  in  terms  of  the  intensity.  /.  and 
phase.  <p.  as  A  -  I'  2(r,z)exp(i<p(r.z)) .  When  the  laser  beam  propagates  in  a  region 

of  uniform  refractive  index,  i.e..  outside  of  the  wiggler,  M2  is  constant.  Simulations 
using  the  FEL  code  MEDUSA  indicate  that  the  beam  quality  of  an  optically  guided 
FEL  beam  is  fairly  good.  For  parameters  similar  to  those  in  Table  I  and  II .  simulation 
results  indicate  the  M 3  »  1 .5  at  the  exit  of  the  wiggler  [26] 


ix)  Distance  Between  First  Relay  Mirror  and  Wiggler 

To  avoid  damage,  the  distance  between  the  wiggler  and  relay  mirror  must  exceed 

(  y.v  V/2 


J  relay 


u- 


(') 


damage  ) 


sinl/20  , 


(9) 


where  0  is  the  tilt  angle  (Fig.  4),  ZR  =  n  R[  /  X  is  the  Rayleigh  range  at  the  exit  of  the 
wiggler.  {/)  =  2(P)tx  R[  is  the  average  intensity  at  exit  of  the  wiggler.  and 

(/damage)  =  50  kW/cm:  is  the  average  intensity  damage  threshold  of  the  mirror.  Using 

the  parameters  in  Tables  I.  II,  and  III,  we  obtain  {see  Table  IV)  for  an  unpinched  optical 
beam. 


and  for  a  pinched  beam 


For  operation  at  A  -  1  fjm  , 


L 


relay 


J  relay 


^  relay 


15m.  0  =  45“ 

8m.  0  =  5° 

5nt.  0  =  45“ 
3m.  0  =  5“ 


27  m.  0  =  45' 

10  m.  0  =  5°* 


and  for  an  unpinched  beam,  and 


(10) 


(ID 


(12) 


9 


11m,  0  =  45’ 
4m,  0  =  5° 


(13) 


L 


relay 


for  a  pinched  beam  (see  Table  V). 


x)  Frequency  Detuning  and  Efficiency  Enhancement 

In  an  FEL  the  efficiency  can  be  increased  by  operating  at  a  frequency  that  is 
detuned  from  the  resonant  valuer^  -2k  d  X.  In  general,  however,  operating  off- 
resonance  leads  to  a  reduction  in  the  growth  rate  [1],  If  the  operating  frequency  is 
denoted  byry,  the  modification  to  the  present  analysis  is  the  replacement  of  the 
wavenumber  shift  tsk  with  A k  - (co-coR)/(2cy:)  in  the  expression  for  the  efficiency  in 
Eq.  (2)  and  in  Eq.  (A7)  which  leads  to  the  dispersion  relation.  Note  that  the  left  hand  side 
of  Eq  .  (A8)  is  not  changed, 

To  demonstrate  these  effects.  Fig.  6  shows  (a)  the  efficiency,  (b)  the  gain  length 
and  (c)  the  filling  factor  as  functions  of  the  relative  frequency  shift  (ru  -a)K)l  coR  for 
X  =  1  pm .  The  minimum  gain  length  is  obtained  near  the  resonance  frequency,  as 
indicated  in  Fig.  6(b),  The  efficiency,  on  the  other  hand,  is  observed  to  rise 
monotonically  as  the  frequency  is  reduced  relative  to  the  resonant  value.  There  is  also  a 
reduction  in  the  filling  factor  for  increasingly-negative  frequency  shifts.  That  is, 
operating  at  lower  frequencies  reduces  the  guiding  effect  and  thus  leads  to  an  increase  in 
the  optical  spot  size. 


IV.  Discussion  and  Conclusions 

Conceptual  design  parameters  for  a  compact  MW-class  FEL  amplifier  operating 
in  the  optically  guided,  high  gain  regime  are  presented.  The  design  presented  here 
assumes  that  the  FEL  amplifier  is  driven  by  a  high  quality  RE  linac  generated  electron 
beam.  Optical  guiding  allows  the  interaction  length  to  be  sufficiently  long  to  achieve 
substantial  growth  of  the  input  signal.  Electron  pulse  slippage,  which  can  limit  the  FEL 
interaction  length,  is  shown  to  be  significantly  reduced  for  a  high  gain  amplifier.  By 
pinching  the  optical  beam  at  the  wiggler  exit  the  Rayleigh  range  is  shortened,  thus 
reducing  the  required  distance  to  the  first  relay  mirror.  The  optical  beam  can  be  pinched 
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by  focusing  the  electron  beam  with  external  magnetic  fields  or  by  utilizing  the  betatron 
oscillation  of  the  electron  beam  envelope.  Grazing  incidence  geometry  permits  a  further 
reduction  in  intensity  to  avoid  damage  to  the  mirror. 

We  find  that  operation  at  2, 141  pm  is  advantageous  because  of  propagation  and 
generation  considerations  as  well  as  for  eye  safety  reasons.  Fora  I  kA.  81  MV  electron 
beam,  a  wiggler  field  of  5  kG.  and  a  frequency  detuning  of  - 1 .5% .  the  power  gain 
length  is  14  cm.  the  intrinsic  efficiency  is  -  1.2%.  the  optical  output  power  is  1  GW.  and 
the  average  output  power  is  1.5  MW.  assuming  a  duty  factor  of  1 .5x  IQ'3  for  the  RF  linac. 
Focusing  the  electron  beam  radius  down  by  a  factor  -  3  at  the  exit  of  the  -  1 .5  m  wiggler 
reduces  the  optical  beam  Rayleigh  range  by  a  factor  of  -  4.  Depending  on  the  mirror 
grazing  angle  and  the  optical  beam  quality,  the  distance  from  the  wiggler  exit  to  the  first 
relay  mirror  can  be  <  3  m.  The  length  of  the  optical  system,  which  includes  the  amplifier 
section  and  the  first  relay  mirror,  can  be  <  5  m. 

A  recent  experiment  at  the  Brookhaven  National  Laboratory  (BNL)  has 
demonstrated  the  optically  guided,  high-gain  FEL  amplifier  concept.  This  experiment 
utilized  a  10  m  long  wiggler  and  an  S-band  RF  linac  that  produced  a  1 00  MeV.  300A 
electron  pulse  of  radius  0.2  mm.  The  wiggler  period  and  strength  parameter  were  3.9  cm 
and  1.1.  respectively.  Experiments  showed  optical  guiding  in  a  high  gain  FEL  amplifier 
operating  at  k  =  0.8  ju  m  using  a  Ti- Sapphire  laser  as  the  seed.  Dipole  steering  magnets 
located  along  the  wiggler  were  used  to  deflect  the  electron  beam  by  ~  lmrad  and  the  light 
beam  was  shown  to  be  steered  and  guided  by  the  electron  beam.  The  transverse  profile  of 
the  optically  guided  light  beam  was  shown  to  be  less  than  half  the  size  of  an  unguided 
diffracting  light  beam.  The  power  peak  output  power  was  -1 50  MW  [8].  Our  simulations 
of  the  BNL  experiment  predict  a  power  e-folding  length  of  42  cm.  a  guided  spot  size  of 
0.24  mm  (1/e  of  intensity),  and  a  peak  power  output  of -ISO  MW,  in  close  agreement 
with  the  experimental  results. 

Simulations  of  the  FEL  amplifier,  taking  into  account  emittance.  energy  spread, 
and  non-uniformity  of  the  optical  guiding  due  to  variation  of  the  electron  current  within 
the  pulse,  are  required  to  fully  evaluate  the  present  design.  Also,  w:hile  estimates  indicate 
that  a  modest  external  magnetic  field  is  sufficient  to  focus  the  beam,  more  refined 
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analysis  is  required  to  design  the  pinch  region  of  the  wiggler.  In  the  absence  of  optical 
pinching,  the  distance  to  the  first  relay  mirror  can  be  <  10  m. 
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Appendix 


The  source  dependent  expansion  (SDE)  formulation  [23]  can  be  used  to  analyze 
the  dynamics  of  the  optical  beam  in  FELs.  A  summary  of  the  principal  SDE  equations  is 
given  in  this  appendix.  For  a  planar  wiggler  the  electric  field  of  the  radiation  can  be 
taken  to  be  linearly  polarized,  and  given  by 

E(r,z)  A(r,z)  exp(ico{z/c  -  ())  +  c.c.  (Al) 

For  the  fundamental  Gaussian  mode  the  complex  envelope  E:i(r,z)  is  written  as 

A(r.z)  =  A„  exp(/4'(z))exp(-(l  -  ia(z))r2 /Rj\z)),  (A2) 

where  An  is  the  amplitude  on  axis. 


'E(z)  =  j*  (M(z')-  jT(z')  )dz' ,  (A3) 

A£(z)  is  the  wavenumber  shift.  T(z)  is  the  spatial  growth  rate.  a(z )  is  related  to  the 
curvature  of  the  optical  wavefronts  and  RL{z )  is  the  spot  size  of  the  beam.  To  allow  for 
an  electron  beam  with  a  radius  Rb  (z)  that  varies  along  the  propagation  axis,  the  electron 
beam  density  is  taken  to  be 

»„(z.r)  =  nh {0)(R2h  (0) / R2h  (r))exp(-  r2  /  R2h  (r)).  (A4) 

The  SDE  equations  for  the  spot  size  and  a(z )  are  given  by  [3] 
dRL(z)  2  ccc(z)  1 


dz 


CO 


R,(z) 


=  -H,  (z)R,  (z), 


and 


dz 

respectively.  Here, 


v  c 

H{z)  =  4—— 


CO 


2  Ir2  ^ 


K1  k 


y  o  l  1  +  K 1 12 


f(z) 


RUZ)  0  +  2  f(z))2  (A k{z)  ~  /T(z))2 


(A5) 


(A6) 


(A7) 
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where  v  =  Ih [kA]/17  is  Budker’s  parameter,  kw  =  2tt/Av,  f{z)  =  R2b{z)f  R\{z)  is 

the  filling  factor  and  the  subscripts  R  and  I  denote  real  and  imaginary  parts.  The  final 
SDE  equation  for  the  phase  shift  and  growth  rate  is  given  by 

♦ Y-T  - 

\0)j  Rj  (r)  R,(z )  dz  2  dz 


(A8) 

Equations  (A5),  (A6)  and  (A8)  describe  the  optical  beam  in  the  high  gain  FEL  amplifier. 
The  electron  beam  can  be  focused  to  produce  a  pinched  optical  beam. 
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■  Aerosol  scattering  [knr1] 

■  Aerosol  absorption  [0. 1  knr1] 

■  Water  vapor  absorption  [0.1  knr1] 


1.625  2.141 


Wavelength  [jim] 


Figure  1:  Aerosol  scattering,  aerosol  absorption,  and  water  vapor  (molecular)  absorption 
coefficients  used  for  the  model  maritime  environment  of  the  propagation  simulation. 
Aerosol  scattering  and  absorption  are  calculated  using  ANAM,  water  vapor  absorption  is 
calculated  using  MODTRAN4. 
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Aperture 
D0  =  80  cm 


Focused  Laser  Beam 
wind 


Target  area 
1 00  cm2 

l::=H 

(p 

\  target 


Range  (focal  length),  L  =  5  km 


Figure  2:  Schematic  of  laser  and  target  configuration  used  in  the  full 
scale  propagation  simulations. 
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Power  on  Target,  (Pm&el)  [MW] 


Figure  3:  Average  power  on  target.  (Patgc[ y ,  versus  transmitted  power,  P ,  in  a 
maritime  environment  for  the  wavelengths  X  -  1.045,  1.625,  2.141  pm  .  Initial  beam 
profile  has  R 0  =  50  cm.  D0  =  80  cm  .  Simulation  geometry  is  shown  in  Fig.  2.  Target 
range,  L  =  5  km;  beam  focus  =  5  km;  target  area  =  100  cm2;  wind  speed,  Vw  =  5 
m/sec;  turbulence  strength,  C2  =  10"15  m'2/3 ;  pointing  jitter  angular  amplitude  =  2 
prad  (white  noise). 
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L 


relay 


Figure  4.  Schematic  of  high  gain  FEL  amplifier  with  a  grazing  relay  mirror.  The  input 
signal  can  be  obtained  a  low  average  power  FEL  oscillator.  The  radiation  beam  is 
optically  guided  in  the  wiggler  and  optically  pinched  at  the  exit.  The  pinched  optical 
beam  has  a  shortened  Rayleigh  range  and  undergoes  rapid  diffraction  upon  exiting  the 
wiggler.  Employing  a  grazing  incidence  configuration  the  resultant  footprint  on  the  relay 
mirror  can  be  made  sufficiently  large  to  avoid  damage. 
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Figure  5.  Plots  of  optical  beam  radius  ( R, )  and  electron  beam  radius  ( Rh )  versus 
propagation  distance  z  .  Wiggler  entrance  is  at  z  =  0  and  the  exit  is  at  z  =  Lw  ~  1.5m. 
The  FEL  operates  in  the  high-gain  regime,  where  the  optical  beam  maintains  a  constant 
spot  size  by  optical  guiding.  Outside  the  wiggler  the  optical  beam  diffracts  and  spreads 
transversely,  as  indicated  in  (a).  By  focusing  the  electron  beam  at  the  wiggler  exit  the 
optical  beam  can  be  pinched  down  to  jndl  =  0.2 1  mm  as  indicated  in  (b).  The 

pinched  optical  beam  diffracts  more  rapidly  compared  to  the  unpinched  example,  i.e.,  has 
a  shorter  Rayleigh  range,  ZR  pinch  =  n  Ripmch  IX  =  6.7  cm.  In  generating  this  figure  a 
frequency  detuning  of  Aro/w  =  - 1.5%  was  used. 
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Figure  6:  Plots  of  (a)  efficiency,  (b)  power  gain  length  and  (c)  filling  factor  as  a  function 
of  frequency  detuning.  The  operating  frequency  is  denoted  by  co . 
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Table  1.  Wiggler  Parameters 


Resonant  Wavelength,  X 

1  |im 

2.141  pm 

Wiggler  Field.  Bw 

8.6  kG 

5  kG 

Wiggler  Parameter.  K 

1.6 

1.86 

Period.  Xw 

2  cm 

4  cm 

Length,  Z,w 

1  m 

1.5  m 

Table  II.  Electron  Beam  Parameters 


Resonant  Wavelength,  X 

1  pm 

2.141  pm 

Energy,  Eh 

76  MeV 

81  MeV 

Current.  Ih 

2kA 

1  kA 

Radius,  Rh 

0.025  cm 

0.02  cm 

Normalized  Emittance,  snFm> 

1 5  mm-mrad 

15  mm-mrad 

Duty  Factor,  D 

10_3 

1.5  xlO'3 

RF  Linac  Frequency 

750  MHz 

750  MHz 

Micro  Pulse  Duration 

1 .3  psec 

2  psec 

Betatron  Wavelength, 

XH=2yXJK 

3.7  m 

7  m 
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Table  111.  FEL  Optical  Beam  Parameters 


Operating  Wavelength.  X 

1  pm 

2.141  pm 

2.141  pm 

Frequency  Detuning.  A  cola) 

0 

0 

-1.5% 

Spot  Size  in  Wiggler.  R , 

0.26  mm 

0.33  mm 

0.4  mm 

Filling  Factor,  f 

0.9 

0.36 

0.25 

Rayleigh  Range  (unpinched,  at  wiggler  exit),  ZR 

22  cm 

16  cm 

23  cm 

Gain  Length  (power).  L(, 

9  cm 

14  cm 

14  cm 

Efficiency  (uniform  wiggler).  rj 

0.8% 

0.77% 

1 .2  % 

Avg.  Intensity  (unpinched,  at  wiggler  exit).  {/) 

1.1  GW/cm2 

0.53  GW/cm2 

0.58  GW/cm2 

Pinched  Laser  Spot  Size,  jnch 

0.1 1  mm 

0.17  mm 

0.21  mm 

Rayleigh  Range  (pinched,  at  wiggler  exit).  ZR  pjnch 

3.6  cm 

4  cm 

6.7  cm 

Average  Pinched  Intensity.  {/pini;h ) 

6.7  GW/cm2 

2.2  GW/cm2 

2.1  GW/cm2 

Peak  Output  Power,  P 

1.2  GW 

0.63  GW 

1  GW 

Average  Output  Power,  (P  ) 

1.2  MW 

1  MW 

1.5  MW 

Number  of  Power  e  -  Folds,  Lv  /  Le 

11 

10.6 

10.5 

Average  Input  Power.  {^npil, ) 

18  W 

43  W 

40  W' 

Optical  Beam  Quality.  M : 

1.5 

1.5 

1.5 
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Table  IV.  First  Relay  Mirror  Parameters  (for  X  =  2.141  pm .  Aoi!  tu=  -1.5%) 


Tilt  Angle,  0  [degrees] 

45  s 

5° 

Distance  from  Wiggler  (unpinched),  Zrelay  [m] 

14 

5 

Distance  from  Wiggler  (pinched),  Irelav  [m] 

8 

3 

The  average  intensity  damage  threshold  was  taken  to  be 
(/*-*)  =  50kW/cm2 


Table  V.  First  Relay  Mirror  Parameters  (for  X  =  1  pm,  Aa>/  co=  0) 


Tilt  Angle,  0  [degrees] 

45° 

5 

Distance  from  Wiggler  (unpinched),  LnlK>  [m] 

27 

10 

Distance  from  Wiggler  (pinched),  Lnky  [m] 

11 

4 

The  average  intensity  damage  threshold  was  taken  to  be 
{l*J)  =  50kW/cm2 
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